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Abstract 

We discuss the neutrino oscillation physics of a very long neutrino factory baseline over a 
broad range of lengths (between 6 000 km and 9 000 km), centered on the "magic baseline" 
(~ 7 500 km) where correlations with the leptonic CP phase are suppressed by matter 
effects. Since the magic baseline depends only on the density, we study the impact of 
matter density profile effects and density uncertainties over this range, and the impact 
of detector locations off the optimal baseline. We find that the optimal constant density 
describing the physics over this entire baseline range is about 5% higher than the average 
matter density. This implies that the magic baseline is significantly shorter than previously 
inferred. However, while a single detector optimization requires fine-tuning of the (very 
long) baseline length, its combination with a near detector at a shorter baseline is much 
less sensitive to the far detector location and to uncertainties in the matter density. In 
addition, we point out different applications of this baseline which go beyond its excellent 
correlation and degeneracy resolution potential. We demonstrate that such a long baseline 
assists in the improvement of the #13 precision and in the resolution of the octant degeneracy. 
Moreover, we show that the neutrino data from such a baseline could be used to extract 
the matter density along the profile up to 0.24% at la for large sin 2 26 l 13 , providing a useful 
discriminator between different geophysical models. 
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1 Introduction 



Input from experiments is central to the search for a theory beyond the Standard Model 
of elementary particle physics. Given the extraordinary robustness of this model, such 
input necessarily requires precision measurements of the relevant parameters in order to 
reveal discrepancies, which, in turn, would provide clues to a higher theory. It is evident 
that at the present time, in addition to the upcoming Large Hadron Collider [1], neutrino 
physics provides an unprecedented opportunity to move successfully towards this goal via 
increasingly accurate measurements of neutrino mass squared differences and mixing pa- 
rameters. These parameters have only been partially revealed by experiments so far (see, 
e.g., Refs. [2-5]). A large number of experiments are thus planned or currently under way to 
achieve this goal. Among these, in the near future, are neutrino beam experiments [6-8] and 
reactor experiments [9-12]. What may be learnt from these experiments depends, however, 
on how large the mixing angle 9 13 is. In particular, some of the important issues crucial 
to the formulation of a unified theory, such as the value of the CP violating parameter 
Sep, or the nature of the neutrino mass hierarchy, may not be determined by these exper- 
iments. With present bounds (sin 2 #13 = Olg 047 , Sqp — [5]), already restricting #13 
fairly stringently, it could, in principle, be very small. This would inevitably result in the 
determination of important physics being relegated to an advanced future facility, such as 
a neutrino factory [13-15]. 

The unprecedented reach and accuracy of a neutrino factory, and the broad scope of physics 
that can be explored has been discussed in detail in several studies (see, e.g., Refs. [16-23]). 
Much of the science centers around isolating the matter effects in neutrino oscillations and 
determining how the various degeneracies and correlations in the measured quantities can 
be successfully resolved to obtain unambiguous physics. The presence of intertwined param- 
eters and effects thus requires careful phenomenological and optimization studies in order to 
best determine the specifications of the source facility and the location of detectors. One of 
the significant results of these efforts is the identification of the "magic baseline" for neutrino 
passage through the Earth, which is a baseline between about 7 000 and 7 500 km [24, 25] . 
Clean measurements of sin 2 26>i 3 and the mass hierarchy are possible due to the disappear- 
ance of CP effects at this baseline. Most recently, these issues have been studied in detail 
from the point of view of the optimization with respect to sin 2 2# 13 , the mass hierarchy, 
and CP violation sensitivities, as well as the measurement of the leading atmospheric pa- 
rameters in Ref. [23]. The conclusion from this study is that an optimal advanced neutrino 
factory would be operated with two detector locations, one in the range between 3 000 km to 
5 000 km, and the second at or around the magic baseline. We recall that the magic baseline 
depends only on the matter density, arising from the condition \^2Gpn e L = 2n, where Gf 
is the weak coupling constant and n e is the electron number density. Obviously, the exact 
density profile plays an important role for the physics potential of the very long baseline, 
whereas the oscillation parameter values and neutrino energies are (almost) irrelevant. In 
comparison, the optimal "short" baseline length depends on the oscillation parameters, but 
the baseline window is large enough such that one does not have to worry too much about 
the actual parameters and the matter density profile for the choice of the detector location. 

Since the precise location of a detector is a function of many variables, a good fraction of 
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which may be unrelated to physics, it is necessary to study the impact of locating the second 
detector not at the optimal magic length, but somewhere in the range between 6 000 km 
and 9 000 km. In addition to this geographical uncertainty, the effects of the matter density 
profile and its uncertainties on the optimization are significant. In this study, we explore 
these questions by introducing a realistic matter density profile model in Sec. 4, and by 
examining the physics consequences of this model in Sec. 5. Note that in our scenario, the 
neutrino factory feeds two detectors, one at L\ = 4 000 km [23], the other at L 2 between 
6 000 km and 9 000 km. We study, as new applications, the 9 13 precision measurement in 
Sec. 6 and the octant degeneracy resolution in Sec. 7, and we sketch the physics case for 
such a baseline in Sec. 9. We believe that the investigation of these questions is especially 
germane in the light of the presently ongoing International Scoping Study for a future 
neutrino factory and superbeam facility [26]. As an interesting application for geophysics, 
we also demonstrate how neutrino measurements may help to extract information about the 
matter density profile and to discriminate between different Earth density models in Sec. 8. 



2 Physics at the Magic Baseline 



It was noticed in Ref. [24,25] that the condition \[2G F n e L = 2ir leads to a disappearance 
of CP violation effects and related degeneracies in P ejX . Several studies have subsequently 
explored the physics possibilities resulting from the enhanced matter effects at and around 
this baseline, see, e.g., Refs. [23,27-32]. A recent closer examination of the phenomenology 
has been done in Ref. [33]. In this section we review the salient features and reasons which 
make this baseline a phenomenologically attractive one. We start with an approximate 
analytical expression for P efM , which is accurate up to second order in the combination of 

two parameters which can (usually) be treated as small, namely a = A ™f ~ ±0.03 and 
sin20 13 [19,29,34]: 

2 „ 2 „ sin 2 [(l-A)A] 
P eM ~ sin 2 2^ sin 2 # 23 - 
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The defining condition for the magic baseline is 

sin(iA) = C ^L = 1X (2 ) 

v2 

for the first (shortest baseline) non-trivial solution. It is apparent that in this case only 
the first term in Eq. (1) remains, which means that any dependence on <5cp is removed. It 
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is important to note that Eq. (2) is also the defining equation for the refraction length of 
a medium [35]. For large matter densities and high energies beyond the resonance energy 
of the medium, the effective neutrino oscillation length in matter approaches the refraction 
length. Moreover, at these baselines and energies, the oscillation length in matter is driven 
by Am^, and not Ara^. When the matter oscillation length equals the refraction length, 
the phase in matter for oscillations driven by the solar mass difference becomes 2-7T, and 
terms containing Am^ and Sqp in the above expansion drop out of the probability [33]. 

We stress that these statements related to the magic baseline are subject to corrections from 
the following: a) The accuracy and the validity of the above expansion in a and sin20 13 , b) 
the validity of the constant density approximation, c) our knowledge of the Earth's density 
profile, and d) the correctness of choosing the average density as the constant density which 
determines the magic baseline via Eq. (2) above. 

The accuracy and validity of various analytical series expansions for matter probabilities 
has been discussed in Ref. [29]. In general, one can choose to expand either in a alone, or in 
sin 26 13 alone, retaining terms up to first order in one of these parameters while treating the 
other parameter exactly. Alternatively, one can choose to treat both as small parameters 
(i.e., do a double expansion), and keep terms up to second order, as we have done above 
in Eq. (1). The single expansion in a retains the exact dependence in sin 29 13, and hence 
is more accurate than the double expansion for values of sin 2#i3 close to the upper bound. 
It is valid for a A = ^j; 1 L <C 1 (which translates to L/E <C 10 4 km/GeV), i.e., when the 
vacuum oscillation length defined by the "solar" mass squared difference is much larger than 
L. On the other hand, the single expansion in sin 2#i3 is most accurate when this parameter 
assumes values which are very small (~ 10~ 3 or less). The double expansion in Eq. (1) is, 
in general, robust over a wide range of parameters (in particular, sin 29\s, L, and E) in the 
sense that relative errors in P^ e , P MT , and P w are restricted to be below 5%. In terms of 

validity, errors resulting from its use are lowest (below 1%) when a A = A ^ L <C 1, i.e., 
L/E <C 10 4 km/GeV, and when sin2#i 3 is small but still above the range which makes 
the single expansion in it preferable. As far as the defining equation for the magic baseline 
Eq. (2) is concerned, note that it is directly obtained from the double expansion. Therefore, 
these constraints on the validity of the expansion need to be kept in mind when drawing 
conclusions about results from a detector located at L — L Magic , especially when precision 
of a few percent is important. 1 

As a first approximation, Eq. (2) can be usefully rewritten as 

W-^-TSOoi^km (3) 
G F n e p(L) 

where 

L 

P(L) = \J P(x)dx (4) 



*An approach independent of an analytical expansion would numerically search for the minimal depen- 
dence on the solar amplitude contributions to the transition probability, but it would still be subject to 
uncertainties in our knowledge of the density profile. 
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is the baseline-averaged density. Solving this equation, we find I/Magic — 7630 km and 
p(7630km) ~ 4.29 gcm~ 32 Most studies in the literature have used p(L) in calculations, 
and have drawn conclusions based on it as the preferred choice for the constant density. 
There is no doubt that the simplified analysis based on the constant density approxima- 
tion has provided many insights into the effect of matter oscillations on neutrinos as they 
travel through the Earth. However, as precision assumes increasing importance in the de- 
termination of neutrino parameters, it is important to ask how well p(L) reproduces the 
profile effect for a given baseline. In particular, we will discuss if p(L) is indeed the best 
choice for the constant density to reproduce physics. We will introduce a constant reference 
density for that purpose which can be different from p(L), and examine its dependence on 
the baseline, on the neutrino oscillation parameters (such as sin 2 2# 13 and 5qp), and on the 
oscillation channel that one is observing. In the next section, we describe the simulation 
methods employed by us, prior to addressing this question in detail in Sec. 4, where we also 
describe our procedure for modeling the Earth's density profile. 

Finally, we stress two important features of the magic baseline relevant to measurements 
of the Earth's matter density. First, the clean dependence of physics at the magic baseline 
on the matter density due to the independence of CP and solar mass-squared difference 
terms suggests that it could be an optimal baseline for Earth density measurements. In 
other words, if the density was considered as just another parameter to be determined 
via neutrino experiments, this determination would be cleanest and most uncluttered by 
correlations if it was done at the magic baseline. Secondly, we note that the condition 
sin(AA) = automatically tells us that the sensitivity to any change in the assumed 
constant density (i.e., the derivative) is maximal here for the terms proportional to sin(AA), 
and consequently, is expected to be high for the full probability. This fact can be usefully 
employed to extract information about the density. These ideas are explored further in 
Sec. 8 of this work. 

3 Simulation methods 

Our simulation is based upon the standard neutrino factory in Refs. [22,23] with a muon 
energy = 50 GeV and a 50 kt magnetized iron detector. We use a total luminosity 
of 4.24 x 10 21 useful muon decays and 4.24 x 10 21 useful anti-muon decays, which can be 
achieved by four years of operation in each polarity at 1.06 x 10 21 useful parent decays/year, 
or by eight years of simultaneous operation of both polarities at 0.53 x 10 21 useful parent 
decays/year/polarity. The neutrino factory uses both the muon neutrino/anti- neutrino ap- 
pearance iy e — > Vy) and disappearance iy^ — > v^) channels. As in Ref. [23], we use a data 
sample without charge identification for the disappearance channels to avoid cuts from the 
charge identification, and therefore increase statistics. For a more detailed description, see 
Refs. [22,23]. 

In order to perform the simulation, we use the GLoBES software [36]. Since it is relevant 
for this study, let us focus on the way GLoBES treats the matter density. In GLoBES, the 

2 Note the non-trivial inter-dependance inherent in seeking a solution to Eq. (3), since p(L) in the equation 
depends actually on ^Magic- 
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matter density of each experiment is implemented as yet another oscillation parameter, i.e., 
the total systematics x 2 f° r one experiment is a function of the oscillation parameters and 
the matter density. For two experiments, the systematics x 2 ys * s then computed as 

Xsys(#i2, 013, 023, $cp, Am^, Aro^, pi, p 2 ) = x!i (^12, 013, 023, ^cp, Aro^, Am 2 31 , pi) + 

+ Xl 2 (0i2,0i3,023,£cp, Am 2 21 , Am^ps) (5) 

for two different baselines Li and L 2 with different density scaling factors pi and p 2 , where 
the oscillation probabilities in x\ x an d xi 2 depend on the respective density scaling factors. 
GLoBES allows arbitrary density profiles, which are scaled by these density scaling factors 
p. For instance, for a constant density profile with density po, the actual density will be 
computed as p ■ po- For an arbitrary profile, p acts as an overall normalization factor with 
which the density in each layer is multiplied. As the next step in the x 2 calculation, the 
external input is added as 

X 2 = X s 2 ys (0i2, 0i 3 , 023, 5cp, A<, Ami, Pi, h) + ( EL ^) 2 + ( EL ^) 2 + • • • ' (6) 

where the dots correspond to other potential external constraints (such as for the solar 
parameters). Typically, we use a p = 0.05, corresponding to a 5% overall normalization 
uncertainty of the density profile. This x 2 is then marginalized over the unwanted oscillation 
parameters and density scaling factors. The marginalization procedure corresponds to the 
projection of the eight-dimensional fit manifold onto the targeted sub-space. For example, 
for a sin 2 2# 13 precision measurement, all other parameters and density scaling factors are 
marginalized over. For a p 2 precision measurement (corresponding to the relative error 
for a constant density), all oscillation parameters and pi are marginalized over, as well as 
a p2 — > oo. 3 Therefore, there is no a priori difference between the oscillation parameters and 
density scaling factors. 

In some cases, we will even use a more detailed approach with different parts of the density 
profile being treated with different scaling factors. In this case, one has to perform the 
density marginalizations in GLoBES manually, i.e., set the density profile in each step and 
manually scan the density scaling factors of the different sub-profiles. For example, we will 
use this approach to measure the lower mantle density only, while we impose an uncertainty 
on the upper mantle density of the Earth. 

For the oscillation parameters, we use sin 2 2# 12 = 0.83, sin 2 2# 23 = 1, Am^ = 8.2 • 
10 _5 eV 2 , Am| x = 2.5 • 10 _3 eV 2 , and a normal hierarchy unless stated otherwise (see, 
e.g., Refs. [2,5,37-39]). In addition, we assume a 5% external measurement for Am^ and 
6 12 from solar experiments (see, e.g., Ref. [38]). We furthermore include matter density 
uncertainties of the order of 5% [40,41] uncorrelated between the different baselines - un- 
less we measure the matter density. In principle, we include all parameter correlations and 
discrete degeneracies [20,24,42,43] where applicable. However, the octant degeneracy will 
not be present for maximal mixing. 

3 Note that in this case, we still impose an uncertainty on pi, because pi is, in the most conservative 
case, fully uncorrelated with p\. 



5 



L = 7500 km 



o 
o 
o 

r 3 



<D 

Q 



2 



f 



Density 

PREM profile 

For simulation: Profile 7 

Reference density: p Re f 



Baseline-averaged density: p 




1000 2000 3000 4000 5000 6000 7000 
x [km] 



4.8 
4.6 

CO 

E 
o 

2 4.2 
4 
3.8 
3.6 



Density mapping: Profile 7 -> p Re f 



p Ref (forsin^2e 13 =0.001,<5 C p=0) 



p 

Parameter dependence 



6000 6500 7000 7500 8000 8500 9000 
L [km] 



Figure 1: Mapping from PREM profile via the simulation profile "ProfikV to our constant reference 
density pR, c f. The left plot shows several possible matter density profiles as function of the propagation 
length x for a fixed baseline L — 7 500 km, as well as the depth as function of baseline. The right plot shows 
our constant reference density pR G f corresponding to Profiler m the left plot as function of the baseline 
L (solid curve) - see text for explanation of the mapping. The baseline-averaged density p is shown for 
comparison as well, and the parameter dependence of the mapping on sin 2 29\z and ^cp is illustrated by 
the gray-shaded area. 



4 Modeling the density profile and determining the optimal con- 
stant density 

Most simulations in the literature so far have used the baseline-averaged matter density 
calculated using Eq. (4), where p(x) is calculated along the baseline using the PREM profile. 
Since the magic baseline is mainly determined by the matter density, this description may 
not be accurate enough to find the optimal detector location (see, e.g., Ref. [44] for a 
discussion at the probability level). Therefore, we follow a different approach in this section, 
and seek a constant density that best simulates the PREM profile. We label this density as 
the "reference density", pR e f- 

To begin with, we note that a numerical simulation cannot use the unmodified PREM profile 
because most numerical techniques are based on the re-diagonalization of the Hamiltonian 
in constant matter density (c/., for instance, Ref. [45]). Therefore, the PREM profile has 
to be made accessible to numerical simulations, i.e., it has to be accurately modeled. One 
possibility is to use a Fourier expansion [46], which, however, does not describe the edges 
very accurately and is not easily implementable in GLoBES. Therefore, we use a different 
approach, which is suitable for the baseline window of interest here (6 000 km to 9 000 km). 
We illustrate this model, which we call'Trofiley" (because of seven layers) in Fig. 1 and show 
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Figure 2: x 2 between a simulated Profile7 and a fit constant density profile as a function of the fit density 
p. The different panels correspond to three different values of sin 2 2#i 3 (see plot labels). In each plot, we show 
three different oscillation channel choices (appearance only, disappearance only and all channels combined). 
The reference density sits at the minimum of the % 2 for all channels combined (chosen for sin 2 20i 3 = 0.001), 
and is marked by the vertical lines in all panels. Note that the scales for the vertical axes differ from each 
other in the three panels. In this figure, 8qp = and L = 7 500 km. 

how it compares to the PREM profile. For any baseline between 6 000 km and 9 000 km, 
the profile looks very similar to the one for L = 7500 km shown in the left panel of Fig. 1. 
We compute Profiler for any given L in this range by choosing the density jumps as edges 
between two layers (such as between upper mantle and lower mantle of the Earth). Each 
layer is then simulated by using the average density given by Eq. (4) within that layer. 
Note that although the use of more steps would make the modeling more accurate, any 
realistic simulation seeks to use an optimal number that retains accuracy while keeping the 
computational time within reasonable limits. We have checked that Profile/ reproduces the 
results of a simulation with significantly more steps with sufficient accuracy. 

As mentioned above, we seek an optimal constant density which can be reliably used as a 
reference density in lieu of full numerical simulations. We define this to be the constant 
density which simulates the PREM profile best and denote it as pR C f. In order to relate 
p Rcf to Profile 7 for any given baseline, we use the x 2 from the complete neutrino factory 
simulation described in the last section (including neutrino and antineutrino appearance and 
disappearance). We identify the pR e f which simulates Profile 7 best by requiring that x 2 be 
minimal between a constant fit profile with p Rc f and the simulated Profiler, for a given set of 
fixed oscillation parameters. This means that we minimize the contribution from the leading 
(zeroth) order profile effect to the total x 2 ■ We illustrate this process in Fig. 2, in which 
we show the variation of x 2 with the choice of the constant density p for different values 
of sin 2 2#i3 and for different oscillation channel choices. The reference density is marked by 
the vertical lines as the minimum of the x 2 between the constant density profile and Profile/ 
(obtained for sin 2 2^13 = 0.001, <5cp = 0. In all three panels (each representing a different 
value of sin 2 2^13 ), the appearance channel dominates the determination of p Re f, but for 
large sin 2 29 1 3 , the disappearance channel also contributes significantly. For small sin 2 26 ) i 3 , 
the appearance channel x 2 becomes much more sensitive to deviations from pR C f than the 
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disappearance channel. The reason is that the disappearance channel is almost in the two- 
flavor vacuum regime described by <-> v T oscillations for sin 2 2# 13 — > 0, i.e., independent 
of matter effects. Therefore, inaccuracies in the estimation of the constant reference density 
are much less important for the disappearance channel than for the appearance channel. 

We show our choice for p Ref (L) as function of baseline as the solid curve in Fig. 1 (right), 
which was obtained for the reference values sin 2 2# 13 = 0.001 and Sep = 0. Expectedly, 
this mapping depends on the choice of oscillation parameters, as illustrated in Fig. 1 (right 
panel) by the shaded area obtained for different sets of parameters. We find that the Scp- 
dependence for large values of sin 2 29 13 is very mild, reflecting the dominance of matter 
effects. As shown above, the minimum x 2 between the profiles p Rc f and Profile 7 is rather 
large (x 2 ~ 15 to 35) and the function x 2 (p) is rather steep for the large sin 2 29i 3 case. On 
the other hand, for small values of sin 2 2#i 3 , the dependence on Sep is stronger, but the 
minimum \ 2 between the profiles p Rc f and Profiler is rather small (x 2 ~ 1 to 10), while 
the function x 2 {p) is more shallow. Therefore, we conclude that the parameter dependence 
in Fig. 1, right (shaded area) may actually not be very significant, either because there is 
no significant 5cp-dependence (large sin 2 2# 13 ), or because the statistical contribution from 
deviations is small (small sin 2 2#i 3 ). We will test this hypothesis later in full numerical 
simulations using different profiles. Note that we will simulate any baseline shorter than 
6 000 km as usual with p. 

A noteworthy point from Fig. 1 is that p Re f is consistently higher than p(L) by about 5%. 
This reflects the dominance of high density regions in the lower mantle (c/., left panel of 
Fig. 1). The fact that matter effects at short baselines are relatively suppressed (irrespec- 
tive of the density profile) has been discussed in Ref. [47]. Thus, in our case, the high 
density central region of the profile contributes with a higher weight compared to the low 
density regions. At short baselines, transition amplitudes for neutrinos can be treated per- 
turbatively, and the leading order terms reflect the contributions of the matter independent 
off-diagonal terms of the effective Hamiltonian, whereas matter contributions reside in the 
diagonal terms. None of this information finds its way into the simple averaging that deter- 
mines p(L). Clearly, an optimized constant density for each channel offers superior accuracy 
compared to simple averaging since it incorporates the effects of matter on transition am- 
plitudes in a weighted way. This has significant consequences for determining the length 
of the magic baseline. For instance, we see from Fig. 1 that the p Ref corresponding to 
7600 km is 4.5g/cm 3 rather than 4.2g/cm 3 . Therefore, the magic baseline length becomes 
an estimated ~ 7300 km according to Eq. (3), i.e., several hundred kilometers shorter. It is 
thus important that prior to fixing the location of the far detector, an optimized reference 
density be used as opposed to using p(L). This sensitivity to the density changes if two 
detectors are used, one at L x — 4 000 km, and the other in the range between 6 000 km and 
9 000 km. We study this combination in the next section. 

5 How precisely does one have to tune very long baseline length? 

In Ref. [23] , the combination of a short baseline L ~ 4 000 km optimal for CP violation and 
the "magic baseline" L ~ 7500 km optimal for degeneracy resolution and mass hierarchy 
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determination was found to have an excellent physics potential. The choice of this specific 
L = 7500 km was (in this and earlier studies [25]) based on the sin 2 29 13 sensitivity being 
optimal there due to vanishing 5 C p effects. While it was found that the short baseline 
is rather insensitive to the specific baseline length in a window between about 3 000 and 
5 000 km, the magic baseline optimum was found in a very small window of L which only 
depends on the matter density profile. In this section, we therefore discuss how precisely one 
has to place the detector on the magic baseline, a question very relevant for the selection 
of detector locations, and the extent of the impact from profile effects and inaccurately 
estimated profiles. 

Before we actually discuss the impact of the long baseline choice on the physics potential, 
let us focus on the most important questions: 

1. Detector location: How much sensitivity does one lose if one moves the detector slightly 
off the optimal length? 

2. Unknown matter density: The length of the magic baseline depends on the matter 
density. What happens if the matter density along the baseline has been misjudged 
within current geophysical uncertainties? 

3. Profile effects: How well does a constant density simulate the matter density profile? 
Does that affect the baseline optimization? Is the actual sensitivity better or worse 
compared to using a constant density? 

In order to answer question 1, we will show the sensitivities in a window between 6 000 and 
9 000km in very small steps. For the complete picture, see, e.g., Ref. [23]. For question 2, 
we will show the results for p^ef as well as the most conservative estimates 0.95 p^ef and 
1.05 pR C f from current seismic wave reconstructions. And for question 3, we will compare 
the results for pR C f with the ones from Profiler. In addition, we will show the curves for p 
to estimate how far earlier studies have been off the more realistic calculations. 

We show in Fig. 3 the sin 2 2#i 3 sensitivity as function of baseline. In the left panel, one can 
clearly see a minimum at about 7300 km for Profiler from a single baseline optimization. 
The reference density pR C f simulates the profile in an excellent way. Misjudging the matter 
density, however, can affect the sensitivity quite severely, and a slightly shorter baseline of 
L ~ 7000 km may be safest from this point of view. 4 However, from Fig. 3, right panel, one 
can read off that while the optimum is still at L 2 ~ 7 500 km, a baseline L 2 ~ 7 SOO!^ km 
does not affect the combined result from two baselines significantly. Additionally, the ex- 
act prediction for the matter density profile is irrelevant. The reason is that the very long 
baseline acts as a degeneracy resolver which can lift the intrinsic and mass hierarchy degen- 
eracies far enough irrespective of the exact baseline choice. Note that placing all detector 
mass at L\ = 4 000km is shown for comparison (horizontal lines), and one can clearly see 
the order of magnitude improvement coming from the second baseline. We have checked 
these results for the combination with L\ = 3 000 km, and also for a larger value of Amjj. 
In both cases, none of the results change qualitatively. In addition, we show the impact of 

4 For a shorter baseline L ~ 7 000 km, the loss in sensitivity is much smaller when misjudging the matter 
density (such as actually having a 5% higher density) than for a baseline fine-tuned to 7300 km. 
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Figure 3: The sin 2 29^ sensitivity as function of L including systematics, correlations, and degeneracies 
at the 5er confidence level. The left plot corresponds to a single baseline only, while the right plot shows the 
combination with a shorter baseline L\ — 4 000 km. In both cases, placing all detector mass at L = 4 000 km 
is shown for comparison as horizonal lines (the index "2L" refers to "double luminosity" because the detector 
mass in the right plot is twice as high as in the left plot). The different curves represent our different matter 
density models as represented in the plot legends. 

a smaller detector mass at L 2 in Appendix A. We conclude that for the sin 2 26*i 3 sensitivity, 
the detector location can be chosen quite freely if there is sufficient statistics from a shorter 
baseline. 5 

Assuming the sin 2 2# 13 sensitivity as our primary criterion for optimization, one can cross- 
check the results with the mass hierarchy and CP violation sensitivity. Note that an ap- 
proach that assumes mass hierarchy and CP violation as primary optimization criteria would 
be much more complicated because one can optimize for two different degrees of freedom: 
sin 2 2#i3 reach and Sep values ("Fraction of Sep" for which one can discover the mass hierar- 
chy or CP violation). In this study, we choose to use the fraction of Sep for different physics 
scenarios in sin 2 2# 13 to cross-check the optimization. However, unless sin 2 2# 13 > 0.01, it 
will not be possible to predict which scenario is actually realized in nature before a neutrino 
factory is built. Therefore, these physics scenarios can really only serve as cross-checks to 
predict how well the final setup would perform for certain hypothetical choices of simulated 
values. 

It has been demonstrated in Ref. [23] that a very long baseline is absolutely necessary 
for the mass hierarchy determination if sin 2 2# 13 is small. We have therefore checked the 
mass hierarchy sensitivity with this in mind and have found that hierarchy determination 

5 This has implications for the proposed INO detector [48], the location of which will be at distances of 
6560 km from JHF and 7145 km from CERN, thus making it a suitable very long baseline detector. 
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Figure 4: The "fraction of (true) ^cp" for which CP violation can be detected as function of L at the 
3cr confidence level. This number represents the fraction of all possible values of <5cp for which the CP 
conserving values #cp = and tt can be excluded. The different panels correspond to different simulated 
values of sin 2 2#i3, as given in the captions. In all cases, the results are shown in combination with a second 
baseline L\ = 4 000 km, and the reference curves for placing all detector mass at L\ are shows as well as 
horizontal lines. The different curves represent our different matter density models as represented in the 
plot legends. For this figure, a normal mass hierarchy has been assumed. 

should be easily possible for L\ — 4 000 km combined with L 2 in the window between 6 000 
and 9 000 km for any value of <5cp and sin 2 2^13 > 10 -4 at the 3cr confidence level (for 
any of the tested density models). Only at sin 2 26> 13 ~ 10~ 4 does the baseline window 
L 2 = (7300 ± 200) km represent a local optimum (Fraction of 5cp = 100%), but the loss in 
the fraction of Sqp if one moves away from this optimal value is at most 20%. 

It is well known that the short baseline plays a major role in enhancing the sensitivity to 
CP violation. Thus, in Fig. 4, we show the combination of L\ = 4 000 km and L 2 on the 
horizontal axes for the sensitivity to CP violation, where the long baseline serve primarily as 
a degeneracy resolver (note the scale on the vertical axes). For large and medium sin 2 2# 13 
(left and middle panels), one can read off that slightly longer baselines compared to the 
magic baseline are preferred. For example, for large sin 2 26 l 13 , the optimum would be L 2 = 
(7 700 ±500) km. The reason for this pull to longer baselines is the Sqp sensitivity returning 
for baselines longer than the magic baseline, i.e., there is additional statistics on the CP 
violation measurement itself. This preference of longer baselines becomes more dominant 
for smaller values of sin 2 2# 13 , and is also accompanied by a narrowing of the peak. Only 
for very small values of sin 2 2# 13 (right panel), does the magic baseline peak at ~ 7600 km 
become less pronounced because the measurement is statistics dominated and any baseline 
with some CP violation sensitivity helps. Again, placing all detector mass at 4 000 km is 
shown for comparison by the horizontal lines. The increase of the physics performance 
by using the very long baseline is actually best for large sin 2 2(9 13 , where the correlation 
with the matter effect affects the short baseline. Coming back to our main questions, 
Fig. 4 illustrates once more that pR e f models the profile effect very well with respect to the 
positions of the peaks, while the peaks for p are slightly shifted. Note, however, that the 
absolute performance using the matter density profile is up to 1% better in the fraction 
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Figure 5: The (projected) A% 2 as function of log 10 (sin 2 20 13 ) for the true values sin 2 20 13 = 0.01 and 
Sep — 5.89 for different baselines as given in the plot legend. 

of Sqp (profile effect). Misjudging the matter density is only relevant if it turns out to be 
much higher than anticipated and the chosen baseline is very long (effect up to several per 
cent). Taking into account all the information from this section, we conclude that the final 
choice of the very long baseline length may well be determined by the availability of detector 
locations and the storage ring design rather than from physics, because the combination of 
two baselines is rather insensitive to the very long baseline length and the specifics of the 
matter density profile. In particular, a baseline somewhat longer than the magic baseline 
does not downgrade the physics potential. 

6 Application: sin 2$i3 precision measurement 

While it is obvious from the phenomenological point of view that a very long baseline is a 
robust correlation and degeneracy resolver, one could conceivably resolve these degenerate 
solutions by different approaches, such as increased statistics, an improved detector, or 
different oscillation channels. However, in this section, we will discuss an obvious generic 
application of a very long baseline: the precision of 6*13. In general, this measurement 
depends strongly on the true value of 5cp [22]. Since, to leading order, there is no Scp- 
dependence at the magic baseline, it provides a low risk option for detector placement, in 
addition to the possibility of improved precision in 13 . 

We define the sin 2 26* 13 precision as the full width relative error on log 10 (sin 2 26> 13 ), i.e., 



Rel. error on log 10 (sin 2 26>i 3 ) 



log 10 (sin 2 2fl 13 )| 



upper 



log 10 (sin 2 26>i 3 )|i owcr 



(7) 



log 10 (sin 2 20 13 )| 



true 
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where "upper" and "lower" refer to the most upper and lower intersections of the fitted 
Ax 2 with the line A% 2 = 9. We do not include the sgn(Am| 1 )-degeneracy because it is 
hard to define the sin 2 20 13 precision including this information. The chosen value A% 2 = 9 
corresponds to the 3a error for Gaussian errors, which is, strictly speaking, not always 
given in this case. We illustrate our definition in Fig. 5, where the (projected) A% 2 is 
shown as function of log 10 (sin 2 2^13) for a specific set of true values. The error is obtained 
from the intersections with the horizontal line (see arrows). As one can easily see in this 
figure, the (Sep, 6*i 3 )-degeneracy [20] is present in some cases under the chosen A% 2 - It is 
included in the error by our definition. One can also see from this figure that while at the 
magic baseline the (Sep, 6*i3)-degeneracy is non-existent (no <5cp-dependence), it reappears 
for longer baselines (L = 7800 km), but is quickly lifted over the chosen confidence level 
(L — 8 100km). Therefore, we expect jumps in the #13 precision as a function of baseline 
wherever this lifting occurs. 

We show in Fig. 6 the precision of sin 2 2# 13 as function of the very long baseline length 
L. The left column corresponds to the very long baseline only, the right column to the 
combination with a shorter baseline at 4 000 km. The different rows correspond to different 
(true) values of sin 2 26*13 as given in the plot captions. The dependence on the true value 
of Sqp is shown by the bands: The upper ends correspond to the worst case Sep, the lower 
ends to the best case Sep, and the thick curves to the "typical value of Sep" , the median 
of the distribution (the precision will be better in 50% of all cases, and worse in the other 
50%). The gray horizontal bands represent all detector mass at the short baseline. 

Let us first of all focus on the left column of Fig. 6 and discuss the 9 13 precision for a 
single very long baseline. There are two aspects which can be inferred from this figure: 
The absolute performance and the risk minimization with respect to Sep- For the absolute 
performance, compare, for instance, the thick dark curve for the "typical" Sep with the 
reference thick horizontal line for L = 4 000 km. One can easily read off that in all cases of 
sin 2 20 13 (rows) a baseline of around 6 700 km to 7 700 km performs significantly better than 
the short baseline for the typical Sep- For large or small sin 2 2#i3 , this baseline window is 
even larger. As far as risk minimization is concerned, the worst case performance (upper 
ends of bands) is almost always better than for the short baseline, whereas the best case 
performance is only better for large sin 2 26* 13 . The reason is the importance of the matter 
density uncertainties for large sin 2 20 13 [41], which can be reduced by a clean determination 
of sin 2 2# 13 at the very long baseline [23]. In addition, one can clearly see that the dependence 
on Sep is minimal at the magic baseline, where the bands become very narrow. The jumps 
in the best case precision come from the lifting of the (Sep, 6*i3)-degeneracy as discussed 
above. In this case, the precision is determined by a different value of Sep, leading to the 
jump. 

As for the combination of two baselines, we compare the very long baseline combined with 
4 000km with all detector mass located at 4 000km (gray shaded region). In all cases, 
the best case, the worst case, and the median performance are better than for the short 
baseline only, which means that the combination of the two baselines is very synergistic. In 
particular, the performance for the typical Sep is much better for the combination than for 
one baseline only. In addition, there is much less sensitivity to the exact value of the very 
long baseline length compared to one baseline only. Note that the behavior as function of the 



13 






Figure 6: The precision of sin 2 201$ (as defined in the main text) as function of the very long baseline 
length L. The left column corresponds to the very long baseline only, the right column to the combination 
with a shorter baseline at 4 000 km. The different rows correspond to different (true) values of sin 2 2#i3 as 
given in the plot captions. The dependence on the true value of <5cp is shown by the bands: The upper ends 
correspond to the worst case <5cp, the lower ends to the best case Sep, and the thick curves to the "typical 
value of 5cp", ^- e -j the median of the distribution (the precision will be better in 50% of all cases, and 
worse in the other 50%). The gray horizontal bands represent the short baseline only for reference (worst 
case, median, and best case). Note that in the right column, double luminosity was used for the reference 
(corresponding to all detector mass at 4 000 km). In this figure, the reference density /9R e f was taken to 
simulate the matter density profile. 
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Figure 7: Regions where the (#23,7r/2 — t3 2 3)-degeneracy can be resolved as function of sin 2 2t3 13 and 
sin 2 #23 for different baseline options as given in the plot labels. The different shadings represent the la, 
2a, and 3<r confidence level sensitive regions (from light to dark, 1 d.o.f.). The different panels represent 
different baseline configurations. Here the true <5cp = and a normal mass hierarchy is assumed. Note that 
this figure does not take into account the mixed degeneracy in sgn(Am 2 1 ) and (#23, t/2 — #23) combined. 



baseline depends very much on sin 2 2^3: For small sin 2 26*13, shorter baselines are preferred 
because of better statistics, whereas for large sin 2 2#i3, longer baselines are preferred since 
matter effects are dominant. 6 



7 Application: Resolving the (0 2 3,7r/2 — 2 3) _ degeneracy 



Another interesting application of a very long neutrino factory baseline is the resolution 
of the (#23,7r/2 — #23)-degeneracy [43]. For example, at the magic baseline, we have from 
Eq. (1) 

* sin 2 2^ sin 2 g 23 Sin2[(1 - f 2 )A] , (8) 

which should be very sensitive to sin 2 #23 and therefore the octant at least for large sin 2 2# 13 . 
However, for small sin 2 20 13 , -P e ^ aglc — at this baseline, which means that there is no octant 
sensitivity. However, in the limit sin 2 2f5 13 — > 



p L=4000k m „ a 2 CQS ^ 23sin 2 2 g l2 Sin2 ^ A ) (9) 



at the short baseline, which means that the octant degeneracy can, in principle, be resolved 
by the cos 2 f^-dependence in the solar term at the shorter baseline. Note that this term 



6 The reason for the preference of longer baselines may actually be non-trivial: At the magic baseline, 
the sin(j4A)-term in Eq. (1) changes sign. Therefore, for longer baselines, the role of r5cp and —Sqp is 
exchanged, i.e., the risk with respect to r5cp is lowered by the combination of the short and long baseline. 
Note that a similar behavior can be found for the silver channel at the short baseline: For maximal mixing, 
the probability is the same as for the golden channel except from r5cp — > —5cp- 
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is suppressed by a 2 ~ 10~ 3 , which means that the octant resolution will be considerably 
worse than for large sin 2 20 13 . For intermediate sin 2 2# 13 , all terms in Eq. (1) will be present, 
which means that the sensitivity to the (# 23 , n/2 — # 23 )-degeneracy will be highly affected by 
correlations with S C p- We therefore expect that the very long baseline could help to resolve 
these correlations. 

We define that the (# 23 ,7r/2 — # 23 )-degeneracy can be considered to have been eliminated 
if, for a given simulated # 23 , no degenerate solution with 823 in the wrong octant (and the 
same mass hierarchy) fits the original solution at the chosen confidence level. In order to 
find the degeneracy, we marginalize over all oscillation parameters in the wrong octant. We 
show in Fig. 7 the sensitivity to exclusion of the (823, n/2 — # 23 )-degeneracy as a function 
of sin 2 2^13 and sin 2 823 for different baseline options. Obviously, for large sin 2 2#i 3 , the 
very long baseline has the better sensitivity because it is less affected by correlations with 
5cp- For small sin 2 2#i 3 , the short baseline can still resolve the degeneracy, but with a 
very poor reach in sin 2 # 23 . For medium sin 2 2# 13 , there is no sensitivity for either baseline 
because the short baseline is affected by correlations, and the long baseline has intrinsically 
no sensitivity. However, the combination of the two baselines does very well in a wide range 
of sin 2 2#i 3 because the long baseline helps to constrain sin 2 2#i 3 very well. Note that we do 
not show the sensitivity to the mixed (octant and sign) degeneracy in this analysis, which 
would impose the additional requirement that the mass hierarchy be resolved (c/., e.g., 
Ref. [49]). In addition, our results are quantitatively comparable with the ones in Ref. [31] 
for the individual baselines. As far as the dependence on <5cp is concerned, we do not expect 
a qualitatively interesting behavior, see Refs. [31,49]. 

It is possible that an octant degeneracy resolution may be obtained earlier than from a 
neutrino factory by atmospheric plus long-baseline data [49,50], reactor plus long-baseline 
data [9,51,52], or long-baseline plus astrophysical data [53, 54]. 7 For large sin 2 26> 13 , a 
neutrino factory would be very competitive to all of these these methods. For small sin 2 2#i 3 , 
a determination from atmospheric (or possibly astrophysical) data combined with long- 
baseline data could come earlier. However, it would require a megaton-size water Cherenkov 
detector and megawatt-class superbeam upgrade. From a practical point of view, it appears 
likely that only one of these combinations (neutrino factory or superbeam upgrade) will be 
eventually realized. 

Finally, we have examined, as a different approach, the efficiency of the silver channel at 
the magic baseline as a degeneracy resolver. At this baseline, we have (see, e.g., Ref. [29]) 

P^ ic ~ sin 2 2fl 13 cos 2 g 23 Sin2[(1 ~ A)A] . (10) 

{l-Af 

Therefore, using Eq. (8), the ratio P efl / P eT = tan 2 # 23 should determine the octant degen- 
eracy without being spoilt by any correlations and degeneracies. We have tested the silver 
channel for the implementation in Refs. [23, 56] with a 10 kt emulsion cloud chamber and 
4 yr of v T appearance, and we have not found a significant contribution to the two baseline 
combination. However, the silver channel does help somewhat for the very long baseline 
alone for intermediate sin 2 26> 13 . The reason for the marginal contribution is the very low 



7 See also Ref. [55] for atmospheric data alone and large #13. 
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Measurement of pR e f 


Measurement of plm 


True value of sin 2 2#i 3 


la 


2a 


3a 


la 


2a 


3a 


sin 2 2^13 = 0.1 


+0.0024 
-0.0024 


+0.0047 
-0.0047 


+0.0071 
-0.0071 


+0.0025 
-0.0025 


+0.0050 
-0.0050 


+0.0075 
-0.0076 


sin 2 20 13 = 0.01 


+0.0065 
-0.0065 


+0.013 
-0.013 


+0.020 
-0.019 


+0.0069 
-0.0068 


+0.014 
-0.014 


+0.021 
-0.020 


sin 2 2(9 13 = 0.001 


+0.021 
-0.019 


+0.043 
-0.038 


+0.067 
-0.056 


+0.021 
-0.020 


+0.043 
-0.040 


+0.068 
-0.059 



Table 1: Relative precision of the measurement of pR C f (constant density along baseline) and plm (lower 
mantle part of the baseline only) at different confidence levels and for different values of sin 2 2#i3 . The same 
parameters as in Fig. 8 are used. Note that for the plm precision in this table, the upper mantle density is 
assumed to be known/fixed. This constraint is studied in Fig. 9. 

event rate in the silver channel at this baseline, i.e., low statistics (about 176 signal events 
for sin 2 20 13 = 0.1). 

8 Application: Matter density measurement 

As mentioned earlier, the sensitivity to matter density changes at the magic baseline is high. 
Moreover, the fact that matter density uncertainties affect the extraction of sin 2 2(9 13 and 
Sep for large sin 2 2(9i 3 is well known (see, e.g., Ref. [41]). This implies that one can extract 
some information on the matter density as well. Such neutrino oscillation tomography using 
a neutrino beam has, for example, been studied in Refs. [57-60]. In this section, we discuss 
the use of neutrino data for potential geophysics applications as a by-product of a very long 
neutrino factory baseline. We stress however that the primary purpose of such a baseline 
remains its use as a degeneracy resolver for neutrino oscillation physics, and the geophysics 
would be a nice addition coming for free. However, if there are different alternative detector 
locations, the geophysics along the neutrino factory baseline might have some impact on 
the choice of the location. 

We follow the treatment of the matter density measurement in Ref. [59] described in Sec. 3, 
i. e., we use the information from both the short (L — 4 000 km) and the very long baseline to 
reduce the impact of correlations. For the short baseline, we use the mean density in Eq. (4) 
and allow for a 5% matter density uncertainty. For the long baseline, we will use several 
models for the density measurement as described below. Note that all of the oscillation 
parameters are marginalized over, i.e., their uncertainties are taken into account. 

As a first model, let us assume that we measure the constant density along the baseline 
modeled by pR C f {cf., Sec. 4). The precision of this measurement is shown in Fig. 8 as 
function of sin 2 2(9 13 (and in Table 1 for specific values of sin 2 2(9 13 ). For large sin 2 2(9 13 , 
where the <5cp-terms in Eq. (1) are small relative to the first term, the measurement is 
very precise and dominated by a combination of resonance peak position and probability 
suppression at high energies, i.e., the characteristic spectral dependence of the matter effect 
(cf., Ref. [60]). For small sin 2 2(9 13 , however, the measurement is dominated by the a 2 -term in 
Eq. (1), which does not have the resonance information and therefore has a lower precision. 
For medium sin 2 2(9 13 , the 5cp-terms act as a background to the density extraction, which 
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Figure 8: Relative precision of the density measurement of pR C f as function of sin 2 29\^ at the la, 2a, 
and 3cr confidence levels (contours from light to dark). For the baseline, we have chosen L = 7 500km in 
combination with L = 4 000 km used to measure the oscillation parameters, as well as we have assumed 
a normal mass hierarchy and 5cp = 0. The dashed curves correspond to not taking into account the 
correlations with the oscillation parameters, which means that the shown combination is close to optimal. 

means that the signal is least clean and the performance is worst there. Note that at the 
magic baseline only the sin 2 26 l i3-term is present by definition. However, the suppression of 
the other terms by sin (A A) = depends on our knowledge of the matter density, which 
means that any deviations in the matter density make these terms return. 

As an alternative model, one can also think about measuring the density along the baseline 
part in a depth d > 670km penetrating the lower mantle of the Earth (c/., Fig. 1, left, 
where the inner density jumps occur). We call this density the lower mantle mean density 
Plm along our baseline, which corresponds to the top density in Profile?. If we assume that 
we know the upper mantle density profile exactly, we obtain the precisions in Table 1, where 
we also compare them to the measurement of pR e f- For large sin 2 29^, we find that one can 
measure the density up to 0.25% at la or 0.75% at 3a. These precisions are very competitive 
to geophysics and could be used to discriminate among different seismic models. 

For a slightly more realistic modeling of the plm measurement, the obtainable precision will 
certainly depend on the local knowledge of the upper mantle (and crust) density close to 
source and detector. This is illustrated by the depth curve in Fig. 1, left, which relatively 
quickly goes into the lower mantle section. Assuming that Eq. (4) describes the density 
measurement to first order, we suspect that the upper and lower mantle densities are highly 
correlated, i.e., for instance, a lower density in the upper mantle can partly be compensated 
by a higher density in the lower mantle. Therefore, we perform a combined fit and discuss 
the precision of plm as function of the error in the upper mantle (and crust) density. Note 
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Figure 9: Left plot: Correlation between lower mantle density /?lm and normalization of the upper 
mantle/crust density /Sum- The arrows indicate the error on plm f° r no constraint on /5um and a very 
precise knowledge on puu- Right plot: Precision of plm as function of relative external error on pum- The 
dashed lines refer to the no constraint limit. Both plots show the la, 2a, and 3<r errors (from light to dark, 
1 d.o.f.) for sin 2 28^ = 0.1. Here we have chosen L = 7 500 km. 

that we only use two parameters in this model: The density plm, which we want to measure, 
and an overall density normalization pum for the profile in the upper mantle, which we will 
marginalize over. 

The result of this analysis can be found in Fig. 9. In the left plot, the correlation between 
Plm and pum * s shown. Because the baseline only runs a short distance in the upper mantle, 
the upper mantle density normalization needs to deviate substantially from unity to cause 
a major effect in plm (c/., scales on axes). The arrows indicate the error on plm for no 
constraint on pum and a very precise knowledge on pum- In the right plot of Fig. 9, we 
impose an external (Gaussian) la-error on pum which represents the precision we believe in 
the upper mantle density close to the chosen source and detector locations. We show the 
precision of plm as function of this relative external error on Pum- As one can easily read 
off this figure, the impact of the upper mantle-lower mantle density correlation will be very 
small as long as one knows the upper mantle density at the chosen location better than 
about 2%, and even for 5%, representing the worst case, the deviation from the best case is 
not very strong. 

We have also tested the dependence of the measurement on the baseline and we have not 
found a significant change in the precision as long as L > 7000 km. However, not sur- 
prisingly, there is a slight optimum at the magic baseline for medium sin 2 26 , 13 , because 
the #cp-dependence is suppressed there. In addition, we have tested the improved neu- 
trino factory from Ref. [23] and have found a considerably higher precision if the muon 
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Fi guxe 10: Different major potential neutrino factory host laboratories, all corresponding detector 
locations at L = 7 200km (dashed curves, corresponding colors), and several chosen baselines to specific 
potential detector locations (red points/end points of curves) projected onto the Earth's surface (solid 
curves). Along these baselines the depth is given in kilometers at several selected points. 

energy is not lowered to 20 GeV. The reason is that the precision part comes from the 
high energy part of the spectrum. Furthermore, we do not expect a quantitatively very 
different result for different values of 5qp or the mass hierarchy in the high precision limit, 
because the <5cp-dependence is of secondary importance for large sin 2 2# 13 , and we have con- 
sidered symmetric operation of the neutrino factory symmetric with respect to neutrinos 
and antineutrinos. For details of the parameter dependencies, see also Ref. [60]. Note that 
compared to Ref. [60], we find slightly higher precisions for large sin 2 26 , 13 , because we use 
the full spectral information (in 43 energy bins) and the knowledge from the short baseline. 

As far as the geophysical interpretation and the selection of detector locations are concerned, 
we show in Fig. 10 different major potential neutrino factory host laboratories and several 
chosen baselines to specific potential detector locations (red points/end points of curves) 
projected onto the Earth's surface (solid curves). Along these baselines the depth is given in 
kilometers at several selected points. Comparing Fig. 10 with seismic wave reconstructions 
(see, e.g., Ref. [61]), one can read off that the obtainable precision on p LM is sufficient for a 
3cr discrimination of different seismic wave reconstructions at specific places, and possibly 
sufficient for the discrimination of different geophysical hypotheses. However, we leave the 
drawing of stronger conclusions based on these results to experts in geophysics. 

9 Physics case for a very long baseline 6 000 km < L < 9 000 km 

In this section, we summarize the results from this and earlier studies to qualitatively 
review the physics case for a very long neutrino factory baseline. We demonstrate that such 
a baseline is extremely versatile over the full allowed range of values of sin 2 29 is, even if the 
specific advantages it has differ for low, medium and high values of this parameter. Note 



20 



that we focus on the baseline window 6 000 km < L < 9 000 km used for this study. 
All allowed values of sin 2 26 1S 

A very long baseline helps for the precision measurements of Am^ and especially 823, 
where, as a general rule of thumb, the longer the baseline, the better [23]. Furthermore, 
as illustrated in Sec. 7, the ability of a shorter baseline to resolve the (823,71/2 — 823)- 
degeneracy is significantly improved for all values of sin 2 2813. In addition, the very long 
baseline improves the sin 2 28i 3 precision measurement in all practical cases when sin 2 2#i 3 
is large enough such that this measurement is relevant. It also minimizes the uncertainty 
stemming from the unknown value of <5cp (c/., Sec. 6). 

Large sin 2 20 13 > 10~ 2 

For large sin 2 2813, the very long baseline helps to reduce the impact of the matter density 
uncertainties by increasing the accuracy in the measurement of sin 2 2#i 3 , which means that 
the CP violation potential of the short baseline becomes significantly enhanced [23]. Note 
that in this case, an improved knowledge of the matter density profile becomes unimportant. 
Therefore, the very long baseline could be one of the key elements in establishing the 
physics case for a neutrino factory for large sin 2 26 l i3. In fact, the very long baseline can 
even measure the matter density along its length to an extremely high precision competitive 
with geophysical techniques, which may provide additional information on the Earth's lower 
mantle (c/., Sec. 8 and Ref. [60]). 

Medium KT 3 < sin 2 20 13 < 10~ 2 

In this scenario, the very long baseline guarantees the mass hierarchy sensitivity for all 
values of 5 C p, and it resolves correlation and degeneracies affecting the CP violation per- 
formance [23]. In addition, the Sqp precision measurement, which strongly depends on the 
true value of Sqp itself for the short baseline only, becomes a low risk endeavor [62]. Fur- 
thermore, there may still be valuable constraints on the matter density profile (c/., Sec. 8 
and Ref. [60]). 

Small lCr 4 < sin 2 20 13 < KT 3 

In this case, the main usefulness of the very long baseline lies in improving the sin 2 2#i3 and 
mass hierarchy sensitivities [23,25]. In particular, the mass hierarchy discovery is guaranteed 
for all values of 5qp, and the CP violation discovery will be possible for at least 50% of all 
values of Sqp [23] . 



21 



Zero" sin 2 20i 3 (< 10 -4 ) 



Even if sin 2 29i 3 is far below the limit of a neutrino factory and Sep becomes unmeasurable, 
there are several interesting measurements which can only be done with a very long baseline . 
For example, the MSW effect in Earth matter could still be verified at a very high confidence 
level because the solar appearance term is large enough [63]. This term would also allow for 
some information on the octant (see above). In addition, the disappearance data could be 
used towards a possible future measurement of the neutrino mass hierarchy [64]. Finally, 
this case may be very interesting from the theoretical point of view since it points towards 
either an exact or softly broken symmetry. 

10 Summary and conclusions 

We have studied the physics and applications of a very long neutrino factory baseline. Our 
work was primarily motivated by the magic baseline [L ~ 7 500 km) where correlations and 
degeneracies naturally disappear. One of the important questions has been the modeling of 
the matter density profile, because the density is crucial for the determination of the magic 
baseline length. We have found that the optimal constant density describing the physics at 
such a baseline is about 5% higher than the average matter density at this baseline. Note 
that this optimal density also describes profile effects with sufficient accuracy. This implies 
that the magic baseline is significantly shorter than one may naively infer, a conclusion that 
has significant implications for decisions on future detector locations. 

Furthermore, we have have re-investigated the baseline optimization for different matter 
profile assumptions. While a single baseline optimization (magic baseline only) is very 
sensitive to errors in the assumed density profile, we have demonstrated that the combination 
with a shorter baseline L ~ 4 000 km allows for a much larger baseline window for detector 
placement. In particular, for CP violation measurements, much longer baselines L ~ 7500 
to 8 500 km could even be preferable, or, at least, do not harm. 

As far as potential applications of a very long baseline are concerned, the primary pur- 
pose remains degeneracy resolution. Beyond that, we have demonstrated other interesting 
applications. For example, we have shown that such a baseline allows for a low risk (with 
respect to 5cp) sm2 2#i3 precision measurement, which is, in the worst well as on the 

average, significantly better than for a short baseline only. In addition, we have illustrated 
the importance of the very long baseline for the octant degeneracy resolution for interme- 
diate 10~ 4 < sin 2 26 l i3 < 10~ 2 . Finally, we have demonstrated that one can also extract the 
matter density along the baseline with a precision of about 0.24% (lo"). If one only wants 
to extract the average lower mantle density of the Earth along the baseline, we still find 
a precision of about 0.4% for a 5% matter density uncertainty in the upper mantle (and 
crust). 

The physics case for a very long neutrino factory baseline has been summarized in Sec. 9. 
Depending on sin 2 26*13, this baseline has different applications and advantages, but it implies 
rich possibilities for all scenarios. Therefore, we conclude that a very long neutrino factory 
baseline would be an extremely versatile tool for a broad spectrum of physics scenarios, and 
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Figure 11: Same as Fig. 3, right, for different detector masses at the very long baseline L 2 as given in 
the plot legend (the detector mass at L\ = 4 000km is fixed to 50 kt). For the matter density profile, the 
reference density pR e f is used. 

strongly warrants inclusion in the planning of a future neutrino factory complex. 
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A How much detector mass is needed at the very long baseline? 

In order to test the impact of detector mass, we show as one example in Fig. 11 the depen- 
dence of the sin 2 2# 13 sensitivity on the very long baseline detector mass for the combination 
with the shorter baseline. One can easily see that mDct > 20 kt at the very long baseline is 
a safe choice, because smaller detector masses lead to a too large loss of events for baselines 
L 2 > 8 500 km. 
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